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LIFT HYSTERESIS AT STALL AS AN UNSTEADY BOUNDARY-LAYER PHENOMENON 1 * 3 


By Franklin KL Moore 


SUMMARY 

Analysis of rotating stall of compressor blade rows requires 
specification of a dynamic lift curve for the airfoil section at or 
near stall, presumably including the effect of lift hysteresis. 
Consideration of the Magnus lift of a rotating cylinder suggests 
performing an unsteady boundary-layer calculation to find the 
movement of the separation points of an airfoil fixed in a stream 
of variable incidence. Then consideration of the shedding of 
vorticity into the wake should yield an estimate of lift increment 
proportional to time rate of change of angle of attack. This 
increment is the amplitude of the hysteresis loop. 

An approximate analysis is carried out according to the 
foregoing ideas for a 6:1 elliptic airfoil at the angle of attack 
for maximum lift. The assumption of small perturbations 
from maximum lift is made, permitting neglect of distributed 
vorticity in the wake. The calculated hysteresis loop is counter- 
clockwise. The computed increment of lift coefficient is quite 
large, indicating appreciable unsteady lift hysteresis for a very 
small reduced frequency of the flow oscillation . It is assumed 
that to the order of this analysis, the wake begins at the separation 
point defined by zero shear. This assumption is questionable 
for unsteady flow. 

Finally, a discussion of the forms of hysteresis loops is 
presented; and, for small reduced frequency of oscillation, it is 
concluded that the concept of a viscous “time lag” is appropriate 
only for harmonic variation of angle of attack with time at 
mean conditions other than maximum lift. 

INTRODUCTION 

The phenomena of “stall flutter” and “rotating stall,” 
which may appear in an axial-flow compressor, both involve 
fluctuations in flow about blades operating near their aero- 
dynamic stall point; that is, at an average flow incidence 
angle near that corresponding to maximum blade lift. 

The analysis of stall flutter has been held back by uncer- 
tainty as to the dependence of airfoil lift and moment on a 
fluctuating incidence angle near stall. In reference 1 there 
was proposed an assumption of the linear aerodynamic force 
and moment relations appropriate to steady flow at a small 
angle of incidence, modified by the further assumption that, 
as the airfoil oscillates, the forces and moments lag behind 
the angular displacement of the airfoil, owing to viscous 
effects. Such a time lag represents an unsteady hysteresis 


which may provide cyclic work to amplify or maintain 
flutter. Perhaps the first experimental study of lift hysteresis 
was that of Farren (ref. 2). Halfman, Johnson, and Haley 
(ref. 3) and Schnittger (ref. 4) have more recently studied 
aerodynamic hysteresis experimentally and have presented 
empirical analyses of their results. 

An analysis of rotating stall (e. g., that of Sears, ref. 5) 
also requires specification of a dynamic lift-incidence relation 
(or the equivalent, as in the study of Emmons, Pearson, and 
Grant, ref. 6, and in Marble’s analysis, ref. 7). Sears has 
adopted Mendelson’s phase-lag hypothesis, and this phase 
angle is an undetermined parameter of his analysis. 

The concept of viscous time lag is not entirety satisfactory, 
however, partly because the phenomenon itself is un- 
explained, but chiefly because the concept obviously cannot 
describe a lift-hysteresis loop which might occur at a nominal 
condition of maximum lift. 

The phenomenon of aerodynamic hysteresis presumably 
depends, at least in part, on the airfoil boundary layer. 
Also, in this study, hysteresis is taken to be a fundamentally 
unsteady phenomenon, not explainable by consideration of 
the steady or quasi-steady boundary layer. 1 In the present 
report, consideration is given to the incompressible flow field 
about a single airfoil fixed in a flow of oscillating incidence, 
under the assumption of an unsteady but nearly quasi- 
steady * laminar boundary layer. (This sort of boundary 
layer is analyzed in ref. 8.) The analysis of this flow field 
is undertaken in order to gain an understanding of the cause 
of lift hysteresis and to describe its form and (crudely) its 
magnitude for a special airfoil at maximum lift. 

The basic quasi-steady flow to be used in the present 
analysis is provided by Howarth’s analysis (ref. 9) of the 
way the laminar boundary layer (and, hence, circulation) 
about an infinite elliptic cylinder depends on angle of attack, 
applied at a condition of maximum lift. In reference 9, 
assumptions are made under which the result becomes 
quantitatively inaccurate, though both the approach and 
the result are qualitatively instructive. The same limita- 
tions affect the present analysis. 

Holding the airfoil fixed while the flow direction oscillates, 
simulates the accepted picture of rotating stall, in which 
successive blade passages stall progressively along a per- 
fectly rigid cascade. The somewhat different case of an 


1 Supersedes NAOA TN 3571, “Lift Hysteresis at Stall as an Unsteady Boundary-Layer Phenomenon,” by Franklin K. Moore. 1956. 

* Airfoils with lift curves that break sharply at stall may show lift hysteresis In steady flow, a phenomen distinct from that under study herein. 

3 In a quasi-steady flow, quantities vary slowly enough so that steady-state results apply at each Instant of time although slight variations are permitted from one Instant to the next. 



882 


REPORT 1291 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


oscillating airfoil in a uniform stream, which is appropriate 
to the stall flutter problem, is not analyzed in this report. 
However, there is an example in the “oscillating airfoil” 
category which illustrates the considerations underlying the 
present study} namely, the rotating circular cylinder in a 
uniform stream. If a circular cylinder is fixed in a uniform 
stream, it, of course, experiences no lift.' Further, if it is 
given an angular displacement, its lift does not change, but 
remains zero. Thus, this degenerate “airfoil” may be said 
to be in a stall condition, at maximum lift, in fact. Now, if 
the circular cylinder is given a constant angular velocity of 
rotation about its axis, then a circulation develops and an 
aerodynamic force (Magnus force) transverse to the flow 
direction is exerted. If the stream velocity is from left to 
right and the rotation is clockwise, then the force is upward 
(lift). If the rotation is counterclockwise, the force is 
downward. 

This phenomenon is explained (ref. 10, par. 27) by con- 
sideration of the boundary layer. In the case of clockwise 
rotation, the upper surface of' the cylinder is moving with, 
and the bottom surface against, the flow. Consequently, 
if circulation remains zero, the boundary layer separates 
later on the top and sooner on the bottom than is the case 
when the cylinder is not rotating. On the top, later separa- 
tion means that the velocity outside the boundary layer is 
lower at separation. Now, the separation point signifies 
the beginning of a wake. Therefore, the clockwise vorticity 
shed into the wake, being proportional to the local outer 
velocity, is less on the top, and the counterclockwise vorticity 
shed at the bottom separation is greater, than in the case of 
no rotation. 

Therefore, owing to clockwise rotation, a net increase of 
counterclockwise vorticity is shed. By the law of con- 
servation of circulation, the circulation therefore cannot be 
zero, and a clockwise circulation must develop about the 
airfoil to compensate for the shed vorticity. According to 
classical hydrodynamics, this circulation results in lift. 

If, instead of rotating steadily, the cylinder undergoes 
a rotational oscillation, the same considerations apply, if the 
reduced frequency of oscillation is small. In that case, the 
oscillating lift is proportional to the instantaneous velocity 
of rotation. Thus, when the “angle of attack” of the cylinder 
is increasing, there is positive lift, and when the angle is 
decreasing, there is negative lift; over a complete cycle, 
the curve of lift against angle of attack would be a loop. 

Therefore, the circular cylinder undergoing a rotational 
oscillation exhibits lift hysteresis, by reason of the response 
of the boundary layer to the movement of the surface. 
In the more complicated problem of a noncircular cylinder, 
or airfoil, similar considerations may be expected to apply. 
Of course, in the airfoil problem contemplated in the present 
study, the acceleration of the flow field may be expected 
to provide an additional component of pressure lift, de- 
rivable from consideration of Kelvin’s impulse. 

PRELIMINARY CONSIDERATIONS 

STATEMENT OF PROBLEM 

Consideration is given to the lift of an isolated airfoil in 
the form of an infini te elliptic cylinder with a semichord 


l and a thickness ratio /3, at a stalling angle of attack a to 
a stream of velocity U (see fig. 1). A full fist of notation 
is provided in the appendix. 

While the airfoil position and the magnitude of U aro 
held fixed, the angle of attack a is permitted to vary with 
time. Such a flow may be constructed by allowing a moving 
source Q to approach an airfoil fixed in an otherwise uniform 



Figure 1. — Notation and coordinate system for elliptic airfoil at 
angle of attack. 



la) 

Sketch (a) Airfoil and source. 


stream, the direction of approach being normal to the stream 
direction, as in sketch (a) : 

As long as h, the instantaneous distance from Q to the airfoil, 
is much greater than l, the airfoil finds itself effectively in 
a uniform stream of incidence 


a=a 


(0)_f. 


Q 

2 irhU 


and of magnitude differing from U only to second order in 
Q/2TrhU. The rate of change of angle of attack is 

a= 2^I7 ^ 


where the dot signifies differentiation with respect to time. 
The foregoing model applies qualitatively to the phenomenon 
of rotating compressor stall, if the moving source Q is taken 
to represent the approach of a flow blockage propagating 
along a cascade. 

The present analysis will be carried out as though a is 
a small constant. Actually, if a is quite small, and lughor 
derivatives such as a are negligibly small, the analysis will 
be correct at each instant using the appropriate instantan- 
eous value of a. (This is the first refinement over the quasi- 
steady assumption which uses instantaneous valuos of a 
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itself; or for oscillatory a, the linear term of a Taylor series 
in reduced frequency.) It is clear that the solution of the 
problem to order a provides a measure of hysteresis: Suppose 
an expression for lift is obtained in the form 


0 L =C^+AaC La +aC Li + . . . 


The first two terms are the quasi-steady contributions. 
The term in a provides that, if a is in the process of in- 
creasing with time, then the lift is higher (ass umin g 0 L{t posi- 
tive) than the quasi-steady value. The converse is true if a 
is decreasing. Thus, if a performs an harmonic oscillation, 
the lift curvo is a loop lying to either side of the quasi-steady 
curve of width proportional to the frequency of oscillation. 

In order that the present analysis bear on the question of 
rotating stall, nominal angle of attack must be selected for 
which the airfoil is in a stalled condition. Maximum lift is 
the most simply described stall condition. Accordingly, 
the nominal angle of attack is chosen as that for which the 
lift is a maximum. This selection is made for two more 
compelling reasons: 

( 1 ) The result will tend to isolate the effect of hysteresis, 
inasmuch as no quasi-steady change in lift results from change 
in a about the maximum lift value. Of course, if lift hy- 
steresis is found under a mean condition of maximum lift, 
then the idea of viscous time lag will thereby be shown to 
bo inappropriate. 

( 2 ) Any other assumption would lead to great theoretical 
complication. The analysis is to be a perturbation of quasi- 
steady flow. If, at the nominal angle of attack, change in a 
resulted in a quasi-steady change in circulation, then, to 
the order of the present analysis, induced wake effects would 
require consideration. 

POTENTIAL FLOW 

Outside the boundary layer of the ellipse, irrotational 
incompressible flow is assumed. At each instant/ the ve- 
locity potential on the surface of the ellipse is (ref. 11 , par. 
71) 

<P=IU (1+/9) cos ( 77 — a)— (1) 

-l7T 

where the surface is defined by 


The foregoing description of the potential velocity dis- 
tribution is made on the assumption that the boundary layer 
is negligibly thin everywhere on the ellipse. This assumption 
is usually quite proper ahead of the separation points. The 
assumption that the wake aft of separation does not import- 
antly affect the potential flow is not proper; certainly, this 
assumption is quantitatively poor at maximum lift, especi- 
ally if the potential flow is used to compute separation point. 
However, the results obtained on the basis of this assumption 
are expected to have qualitative validity. 

BOUNDARY-LAYER ANALYSIS 

QUASI-STEADY BOUNDARY LAYER 

As a basis for subsequent calculation of unsteady effects, 
the quasi-steady laminar boundary layer on the ellipse may be 
approximately determined by the K&rm&n-Pohlhausen inte- 
gral method, as improved by Holstein and Bohlen (see ref. 
12 , ch. XH). The differential equation is 

<1Z_F{k). „ dq_ 

as g as 

subject to the initial condition at the stagnation 
( 2 = 0 ): 

k o =0.0770 

where Z=6 2 U/vl, 6 being the momentum thickness, 
function F (k) is tabulated in reference 12. 

Determination of Z and r at maximum lift. — Given the 
velocity distribution q, the growth of the boundary layer may 
be computed from equations (5), the calculation proceeding 
until both the separation points (s and s, fig. 1 ) are reached, 
for which (k) 2 =(k)i=— 0.1567. 

In the present problem the potential velocity distribution 
has not yet been fully prescribed, since T remains unknown. 
Suppose that, for a given a, a value of r is assumed, and the 
boundary-layer calculation is carried out (ref. 9). Then, at 
the top separation point, clockwise vorticity is shed into the 
wake at the rate 


(5a) 

point 

(5b) 

The 


f 5 j 1 , ~ 

Jo U bn dn ~2 {V?)i 


( 6 ) 


2=1 cos y^fil sin tj ( 2 ) 

Along the surface, measuring s clockwise, 

ds 

^=— (sin 2 77 + 73 * cos 2 77 )“^ = — (3) 

g ^U = R [( 1+/3 ) S1D ( 7 ?~ a )+2 vUl\ 

At this juncture, the condition of maximum lift has not been 
imposed, and the circulation r is left unspecified. Of 
course, both a and T may vary with time. 


while at s, counterclockwise vorticity is shed at the rate 

7 ; («?)»• Now, if +){ differs from («i)j, there is a net change 

of circulation in the wake, which is impossible in steady flow. 
Therefore, new trial values of r must be assumed until the 
particular value of r is determined for which the solution of 
equations (5) yields separation points for which (u,) 2 = («i)i. 

The foregoing procedure must be carried out for each of a 
number of angles of attack in order to determine the values of 
a and T at maximum lift. In reference 9, Howarth per- 
formed these calculations for /3=l/6 and determined a the- 
oretical lift curve. The maximum value of T/2vUl was found 
to be 0.0761 at an angle of attack a=7°. His complete dis- 
tribution of Z at this condition is not presented in reference 9 . 


(The quasl-stoody assumption applies precisely for the calculation of the unsteady velocity potential 
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Accordingly, the calculation has been repeated by integrat- 
ing equation (5a) in the form 


where 


K=a(v)Z 


A(j?) = a(v)=—^ 


and, from equation (4), 

2=^ [| sin 0,-7°)+0.076l] (8) 

At the forward stagnation point where 2=0; 170 = 190.74°, 
and the initial condition is, from equations (5), (6), and (7a), 

Z 0 =0.0770/a o =0.00407 (7c) 

The solution was obtained using a step-by-step method in 
which a parabola is passed through two known values and the 
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Fiqube 2. — Distribution of thickness parameter Z around ellipse. 


Differentiating equation (7a) yields 
< ^=B(rj)Z a +b( V ) 


where 


B=^F'(k) 


bm-A (^F+~pF^) 


The initial condition for Z a at tj 0 is determined by specifying 
that dZajd-q must be finite there; from equations (9a) and 
(9b), 


(7 \ _/ z n<L*-Zq°*F' \ _ 
{Za)0 \ q,F' 


-0.00660 


In equations (9) all quantities are to be evaluated at tho 
condition of maximum lift; the appropriate superscript (0) 
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next unkn own value of dZ/di 7 , integrating to find Z in terms of 
the unknown dZ/di 7 , then applying equation (7a) at the un- 
known point to solve for dZ/dq. Two starting values were 
found from a Taylor series about 170 - The step size in 17 was 10 ° 
except near the stagnation and separation point where finer 
spacing was used. The solution of Z is shown in figure 2 and 
in table I. The separation points, for which *=—0.1567, 
occurred at 17 = 80 . 0 ° and 340.83°. Of course, 2 should be 
the same at s and s. The difference cited in the table indi- 
cates the degree of error present in the calculations. 

Determination of bZfba. — For subsequent use in the 
unsteady equations, it is necessary to know the rate of 
change of Z with a in quasi-steady flow. At maximum lift, 
when a is changed, the quasi-steady boundary layer changes, 
and the locations of the separation points are changed. Of 
course, the velocities at s and s must remain equal, because 
£>T[da=0 at maximum lift, by definition. 


is omitted for brevity. Of course, for purposes of finding 
2 a and 2 a, in equations (9), the angle of attack of 7° should 
be replaced by a in equation (8) and set equal to 7° again, 
subsequent to differentiation. Equations (9) have been 
integrated to yield Z„, by the same method as described for 
finding Z, and the result is shown in figure 3 and in table I. 

UNSTEADY BOUNDARY LAYER 

The next step in the analysis is to determine the depend- 
ence of the boundary layer on the angular velocity a, as- 
sumed small. To this order of approximation, equation (8) 
describing the potential flow must be modified to include the 
possibility of a contribution to circulation in proportion to 
a (or, in dimensionless form, e = al/U), as follows: 

2 [| sin (4-7°)+0.0761+T £ ] (10) 
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The coefficient of circulation hysteresis 7 must be found from 
a condition of vortieity shedding at the separation points of 
the unsteady boundary layer. The contribution to lift 
proportional to a then follows. Determination of the proper 
vorticity-shedding condition will be discussed in a subse- 
quent section. 

The unsteady form of the K&rm&n momentum equation 
(sec ref. 12 ) is 

< u) 


Steady equations (5a) and (7a) are obtained from equation 
( 11 ) by omitting the term on the right side. An iteration 
procedure might then be adopted: The quasi-steady S* can 
be substituted into the right side of equation ( 11 ), and a 
new solution obtained, to include first-order unsteady 
effects. Using the definitions of reference 12 , and writing 
the time derivative in equation ( 11 ) as 


^(M*)=i|;(M*) 


( 12 ) 


there is found, corresponding to equation (5a), 


dZ F(k) 
ds g 




(13) 


The function /dc) is tabulated in reference 12 , and g is given 
by equation (10). Again, for purposes of finding q a and 
5 U,, the angle 7° in equation (10) should temporarily be 
replaced by a. 

Instead of the indicated iteration, in the present study 
the equivalent procedure is adopted of finding the coeffi- 
cients of the expansion 


Z=Z™+AaZ a +eZ c + . . . 


The coefficients Z (0) and Z a have already been found (eqs. (7) 
and (9) and figs. 2 and 3). The derivative Z, remains to 
be found. Differentiating equation (13) and noting that 
K a =q,Z a +q a ,Z, 


Therefore, the definition of momentum thickness 


0 = r ±(l-2£) dy 
J 0 Ut \ 


shows that, if the velocity u in the boundary layer has a 
part proportional to U\ and a part proportional to e, then 
the part of 6 (and hence of Z=8 2 U/vl) which is proportional 
to e must have a simple pole in u x . Actually, of course, the 
quantity 6 is inappropriate for defining a thickness of a 
profile of the type shown in sketch (b), and the appearance 
of a pole in 6 simply indicates this lack of physical significance. 


/ 



The foregoing considerations suggest that a new variable 
W be defined to replace Z„ : 

W=qZ, (15) 

Substituting equation (15) into equation (14) yields 

(1 + F')W+ZF'q a -F 

2flZ [&+§+£ (g.Z„+g M Z)] (16) 


^~F'(k) ^ Z t +F' Z-5 5 - 

ds 2 2 2 

2 f[f + S + 7> z * +2 “ Z) ] (14> 

In equation (14) and hereinafter, evaluation of quantities 
in the steady state at maximum lift is to be understood, and 
the superscript ( 0 ) is omitted for brevity. 

It is impossible to ensure a finite value of dZJds at the 
stagnation point because of the second-order pole 
— 2 /iZgo/g*. The physical reason for this result is the 
fact that, at the stagnation point of a certain instant, the 
boundary-layer velocity profile will not vanish, as in steady 
flow, but rather will respond to the instantaneous acceleration 
c more promptly than the outer potential flow. A profile 
of magnitude e may thus be expected, vanishing in the outer 
stream as well as the wall as shown in sketch (b) : 


In order that dW/d s be finite at the stagnation point, the 
two poles on the right side of equation (16) must cancel, 
yielding the initial condition 

HjSSI.1-- 1 501,1015 (17) 

where the numerical value is obtained using equations (3), 
(7c), and ( 10 ), and the tables of reference 2 . Inasmuch as 
q,=y/R (eq. ( 10 )), the function W may be split into two 
parts, as follows: 

(18) 

so that (changing to 17 as independent variable) equations 
(16) and (17) provide 


^=C(v)X+c(v) 


(19a) 
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CO?) =(1+^0 1 

civ) =2fiRZ (q^ a + qav Z)J 

and 

X 0 = —0.001015 
^=C(v)Y+d(v) 

where 

F R 

d( v )^-ZF'^ 

and 

Fo = 0 


(19b) 


(19c) 

(20a) 

(20b) 

(20c) 


Equations (19b) and (20b) are evaluated using the solutions 
for Z and Z a , the tables of reference 12, and equations (3) 
and (10). 

Solutions of equations (19) and (20), obtained by the 
method used to find Z, are presented in figures 4 and 5 and 
table I. 



Figure 4. — Distribution of parameter X of unsteady boundary layer. 



Coordinate on surface of ellipse, ij, deg 


Figure 5. — Distribution of parameter Y of unsteady boundary layer. 


DETERMINATION OF LIFT 

UNSTEADY BALANCE OF SHED VORTICITY 

In order to determine the unsteady pressure lift, the 
coefficient y (eq. (10)) must be determined. In steady flow 
(see the discussion accompanying eq. (6)), the circulation 
r was obtained by requiring that vorticity be shed in equal 
and opposite amounts at the two separation points. 

In the present unsteady problem, the net rate of vorticity 
appearance in the wake must vanish, not only in the quasi- 
steady approximation, but to order a as well, in view of the 
assumption that the quasi-steady circulation is maximum. 
By the classical theorem concerning constancy of circula- 
tion, any net rate of discharge of vorticity into the wake 
must be balanced by a rate of increase of circulation about 
the body. If the airfoil is no min ally at maximum lift, 
then the circulation terms (eq. (10)), to order e^al/U, are 

0.0761+76 

The rate of increase of this expression is at most of order 
a, and therefore there cannot be any net discharge of vor- 
ticity to order a. 

If the airfoil were not at maximum lift, then the expression 
for circulation would contain a term proportional to Act, 
which would change at the rate a, and would have to bo 
balanced by a net rate of vorticity discharge of order a. 
In turn, this distribution of circulation in the wake would 
induce further modifications of the potential flow. Therefore, 
the assumption of maximum lift permits the neglect of the 
induced effects of distributed circulation in the wake. 

Movement of separation points. — In order to effect a 
balance (to order a) of vorticity shed at the separation 
points, the movements of the separation points must bo 
taken into account. The position of the top separation point 
s may be written 

V={vY,+Aa (g)_+e (§;)_+• • • (2D 


The coefficient (J^)_ is obtained from the quasi-stoady 
solution: 

At separation (*=—0.1567), 

(*V\ 

\Wr Vdic/dij/i 

From equations (7), 

bv~~R ( Z 2-"+ Z <&+72 R s,n 2v ) 


^=— (Zq, a +Z a q,) 
whence, holding k fixed at —0.1567, 




Zg„+Z,g,+;~|| sin 2 V 


), 


(22a) 

(22b) 


(23) 
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where Z and Z a may be obtained from table I. 

The coefficient (dTj/d«)j comes from the unsteady solution: 


/ bn/b<\ 
be)-, Vd/c/dj?/; 


From equations (7), 


K t — j) (Z Zq,,^) 


which, upon substitution of equations (10), (15), and (18), 
becomes 


(X+77) sin 2n (24) 


Equations (22a) and (24) thus provide that 



£1+^1 

Zg„+Z, & +^ sin 2 V 


(25) 


At the bottom stagnation point, equations (21), (23), 
and (25) apply, with subscript 8 replacing subscript a. 

Equation of vorticity shed at separation points. — Taking 
into account the motion of the separation point, the rate of 
vorticity shedding at the top is given by the following 
equation, which replaces equation (6): 



bu 

bn 


j l_ a , d8 - 

dn -~2 Ul l dt Ul 


From equations (3) and (21), 


(26) 



Equations (27) yield the following expression for the right 
side of equation (26) : 


\ (uf)?' +(«i)r T A« 


(^h^v\ 
\b v C to)-. 


+ € 


(tS+ ra S)J <**> 


Expression (28) represents clockwise vorticity shed at the 
top separation point. At the bottom separation point, the 
amount of counterclockwise vorticity shed is also represented 
by equation (28), if subscripts s are replaced by s. 

Therefore, equating the net discharge of vorticity to zero, 



In the quasi-steady flow, (-u 1 )-= («,)„ and the coefficient of 
Aa must be zero. Therefore, the coefficient of e in equation 
(29) must vanish: 


(t S)r (It %\+ u [( b I H* Is).] -0 


or 


/ dg/bT? bij/be \ F / bi]/ba \ / bvjbe \ ~| (R bv/ba), 

V Rbv/bcc AL VbvIbeJAbVjbaJ.J^ 1 (R blj/ba)-~ U 

(30) 


Equations (23), (25), and (30) and table I suffice to deter- 
mine y (which appears in eq. (25)). The result is, 


Y= — 6.1 (31) 

Definition of separation point. — In effect, it has been 
assumed that, dining the unsteady motion, separation is 
defined by the condition of zero shear (k=— 0.1567) and 
the subsequent appearance of reverse flow relative to the 
surface, just as in steady flow. This assumption is open to 
question. The question is how (or whether) local velocity- 
profile characteristics may be interpreted to identify the 
leading edge of a wake. 

The usual steady criterion, which notes the appearance 
of reverse flow downstream of the point of zero shear, implies 
that the fluid in the wake is fixed relative to the body. If, 
in the unsteady case, the wake may still be regarded -as 
fixed to the body, then it may be that the steady criterion is 
still applicable. 

However, the present assumption of the steady definition 
of separation is not advanced with complete confidence. 
Rather, it is felt that only a suitable experiment can settle 
this point. 

LIFT OF AIRFOIL 

The steady lift coefficient of the airfoil of figure 1 is deter- 
mined from the steady circulation: 

CL m = - pU r — =TyT=27r(0-0761)=0.48 (32) 

\pU\2l) 


There are two contributions to lift proportional to e. From 
the unsteady circulation, 

Oz^=27t7= — 38.6 (33) 

and a further contribution is found by consideration of the 
remainder of the potential flow. The two components of 
Kelvin’s impulse for the flow illustrated in figure 1 (leaving 
circulation out of account) are 

R,l, =irpUP(p 2 cos a, sin a) 

(See pars. 71 and 123 of ref. 11.) Whence, the corresponding 
components of vector force are 

F x , F y =a -~^=7rpUl 2 a ( — (3 2 sin a, cos a) 
and the lift is 

L=F cos a—F x sin a=TpUl t a —^r~ ^l+yq^p cos 2a ) 
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The following lift coefficient results : 

<?£=§ (1+JS 2 ) 003 2«)=3-09 (34) 

The final expression for lift combines equations (32), (33), 
and (34): 

Ci=C$>+<:{C2> t +C£l) 

=0.48—36 ^ (35) 

The sign of the second term of the result of equation (35) 
indicates that, while angle of attack is increasing, the lift 
is lower than the quasi-steady value, and higher if the angle 
of attack is decreasing. Thus, near maximum lift, the lift 
curve would exhibit a counterclockwise hysteresis loop 
enclosing the stall point. This result is perhaps counter to 
expectations, because clockwise hysteresis is found experi- 
mentally for oscillating airfoils. It may be that different 
directions of hysteresis should be expected when the airfoil 
oscillates and when, as in the present study, the stream direc- 
tion oscillates. 

In any case, it may be shown that the overriding effect 
producing counterclockwise hysteresis in the present prob- 
lem is the quasi-stead 3 r movement of the separation point 
over the top surface. As this separation point moves forward 
under increasing angle of attack, clockwise vorticity in the 
boundary layer joins the wake as the separation point passes. 
Accordingly, a counterclockwise airfoil circulation (negative 
y) is required to balance this effect. The term of equation 
(26) that is concerned with this movement is the one in- 
volving ds/dt. 

For the elliptic airfoil problem treated herein, the quasi- 
steady movement of the upper separation point is quite 
extensive; numerically, Zhjfda= 13.9, indicating that the 
separation point is very loosely fixed to the airfoil suface. 
By way of contrast, in the case of the rotating circular cy- 
linder, there is no effect of this sort, because a change in 
angle of attack produces no quasi-steady movement of the 
separation point at all. Other contributions to the shedding 
of vorticity then lead to the result of clockwise hysteresis 
for the circular cylinder. 

LIFT HYSTERESIS 

The foregoing analysis does not provide a complete theory 
of unsteady flow about a stalled airfoil. Rather, the analysis 
illustrates the considerations that would underlie such a 
theory and further make plausible the general assumption of 
an expansion of lift coefficient in the form of equation (35). 
This expansion would be valid for nearly quasi-steady condi- 
tions. Also, it has been shown that a counterclockwise 
hysteresis loop may be expected at a nominal condition of 
maximum lift with, apparently, a large amplitude. 

If the foregoing conditions are met, and A a is simple 
harmonic, 

Aa cc sin 2iraf (36a) 

then 


Use of relation (36b) in the lift formula (35) may be said to 
correspond to the assumption (ref. 5) of a positive lift- 
curve slope and a phase lag, 90° in this case, though the 
positive slope would not correspond to the steady lift curve. 
The hysteresis loop for this case appears as an ellipse on the 
lift curve of figure 6 (a). The amplitude and width of the 
loop are assumed small in the present discussion and are 
exaggerated in figure 6. 

If a is not simple harmonic, then the concept of phase lag 
is altogether inappropriate. For example, if a is more 
nearly a “saw-tooth” function of time (fig. 6 (b)), then tho 
lift increment is nearly a “battlement” function. The 
corresponding hysteresis loop is nearly rectangular. If a 
changes according to an exponential pulse (illustrated by a 
.Gaussian curve in fig. 6 (c)), then the hysteresis loop is 
egg-shaped, with the broad end to the right. 




(a) Harmonic a(i). 




(b) 

(b) “Saw-tooth” function a (1). 



(c) Exponential pulse a(<). 

Figure 6. — Variation of lift coefficient and angle of attaok with time 
and with each other. 


a cc— 27i« sin'(2mwf— 90°) 


(36b) 
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CONCLUSIONS 

Tlic analysis of rotating stall in an axial-flow compressor 
requires specification of a dynamic lift curve applicable 
near stall. It has previously been suggested that unsteady 
lift hysteresis is an important characteristic of such a curve. 

Consideration of the familiar experimental fact of Magnus 
lift on a rotating cylinder indicates a theoretical approach 
to the question of aerodynamic hysteresis which, though 
certainly uot definitive, may prove helpful. The accepted 
explanation of Magnus lift is that, if the cylinder is in motion 
toward the left and rotates clockwise, the movement of the 
surface delays boundary-layer separation on the top and 
advances it on the bottom. In steady flow, considerations 
of constancy of wake circulation require that the outer 
velocity at the two separation points be equal. The de- 
layed separation at the top implies a lower velocity (vice 
versa on the bottom), and a compensatory clockwise circu- 
lation must therefore occur. 

The foregoing reasoning is extended to apply to the prob- 
lem of an airfoil of elliptic section in a stream of constant 
velocity but of (slightly) oscillating direction. The airfoil 
is considered to be nominally at maximum lift. This 
assumption, reasonable for unsteady problems at nearly 
stalled conditions, provides an essential simplification. To 
first order in small quantities, the lift (circulation) increment 
duo to the oscillation can depend only on rate of change of 
angle of attack; and, just as in the cylinder case, all induced 
wake effects may be ignored. 


For purposes of computing the amount of vorticit 3 r shod 
into the wake, the separation point is identified as the point 
of vanishing shear, just as in steady flow. It is not clear, 
however, that this assumption is proper. 

Under these various assumptions, the unsteady increment 
in lift coefficient of the ellipse is found to be —36 al/U, of 
which about 92 percent is due to the unsteady movement 
of the separation points and the remainder is due to impulsive 
pressure. This lift increment gives the amplitude of a lift- 
hysteresis loop at maximum lift. The loop is counterclock- 
wise, a result that can be related to the extremely migratory 
tendency of the separation point on the upper surface of the 
ellipse under a change in angle of attack in steady flow. 

Finally, assuming oscillations of low reduced frequency, 
certain observations may be made concerning the shapes of 
hysteresis loops, and the validity of the idea of a viscous 
time lag in connection with unsteady lift: If the angle of 
attack undergoes harmonic oscillation, then the lift incre- 
ment is also harmonic with a 90° phase lead or lag, depending 
on the sign used in the definition of C lt . In this case, the 
hysteresis loop is elliptic. If the angle of attack varies in a 
nonharmonic manner, then the variation of lift does not 
have the same dependence on time, and the idea of time lag 
is inappropriate. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 17, 1955 
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APPENDIX 


SYMBOLS 


The following symbols are used in this report: 

A, a functions of 17 (eqs. (7)) 

B, b functions of 17 (eqs. (9)) 

0 , c functions of 17 (eqs. (19b)) 

Cl lift coefficient 

C L rate of change of lift coefficient with angle of attack a 

C L - rate of change of lift coefficient with a 

d function of 77 (eq. (20b)) 

F universal function for boundary-layer calculation 

(eq. (5a)) 

ji universal function for boundary-layer calculation 

(eq- (13)) 

L lift 

l semichord of elliptic cylinder (fig. 1) 

n coordinate measured normal to surface (eqs. (6), (36)) 
q dimensionless outer velocity, muJU 

R function of 17 (eq. (3)) 

s dimensionless coordinate measured along surface of 

ellipse (fig. 1) 
t time 

U stream velocity (a constant) 

u velocity parallel to surface 

W function related to unsteady boundary layer (eq. 

(15)), -22. 

X function related to unsteady boundary layer (eq. (18)) 
x Cartesian coordinate of surface of ellipse 

T function related to unsteady boundary layer (eq. (18)) 
y Cartesian coordinate of surface of ellipse 

Z function related to momentum thickness of boundary 

layer (eq. (5a)), 

a angle of attack 

a time rate of change of a 

A a increment in angle of attack, =a—a m) 

P thickness ratio of elliptic cylinder 

r circulation in outer flow 

7 coefficient of hysteresis in circulation (eq. (10)) 

5 over-all thickness of boundary layer 

5 * displacement thickness of boundary layer 

€ dimensionless angular velocity, tscd/U 

17 coordinate on surface of ellipse (fig. 1, eq. (2)) 

6 momentum thickness of boundary layer 

k function for boundary-layer calculation (eq. (5a)) 

v kinematic viscosity coefficient 

p density 

r skin-friction coefficient 

<p velocity potential 

oj frequency 

Subscripts: 

s, s evaluation at top or bottom separation point, respec- 
tively, of steady flow at maximum lift (a=a (0) , 
6 = 0 ) 

0 evaluation at forward stagnation point, where 2=0 


1 evaluation at outer edge of boundary layer 
evaluation at lower separation point 
Subscript notation is used for partial differentiation where 
convenient. 

Superscripts: 

(0) steady conditions at maximum lift 

(1) unsteady contribution due to movement of separa- 

tion points 

(2) unsteady contribution due to impulsive pressure 
evaluation at upper separation point 

' denotes ordinary derivatives 

REFERENCES 

1. Mendelson, Alexander: Effect of Aerodynamio Hysteresis on 

Critical Flutter Speed at Stall. NACA RM E8B04, 1948. 
(See also Jour. Aero. Sci., vol. 16, no. 11, Nov. 1949, pp. 
645-654.) 

2. Farren, W. S.: Reaction on a Wing Whose Anglo of Incidonoo is 

Changing Rapidly. R. & M. No. 1648, British A. R. C., 1936. 

3. Halfman, Robert L., Johnson, H. C., and Haley, S. M.: Evaluation 

of High-Angle-of-Attaok Aerodynamio-Derlvatlve Data and 
Stall-Flutter Prediction Techniques. NACA TN 2633, 1961. 

4. Sohnittger, Jan R.: Single Degree of Freedom Flutter of Compressor 

Blades in Separated Flow. Jour. Aero. Sci., vol. 21, no. 1, 
Jan. 1954, pp. 27-36. 

5. Sears, W. R.: A Theory of “Rotating Stall’’ in Axial Flow Com- 

pressors. Graduate School Aero. Eng., Cornell Univ., Ithaca 
(N. Y.). (Contract AF 33(038)-21406.). 

6. Emmons, H. W., Pearson, C. E., and Grant, H. P.: Compressor 

Surge and Stall Propagation. Traus. A.S.M.E., vol. 77, no. 4, 
May 1955, pp. 455-467; discussion, pp. 467—469. 

7. Marble, Frank E.: Propagation of Stall in a Compressor Blade 

Row. Jour. Aero. Sci., vol. 22, no. 8, Aug. 1955, pp. 641-664. 

8. Moore, Franklin IL: Unsteady Laminar Boundary-Layer Flow. 

NACA TN 2471, 1951. 

9. Howarth, L.: The Theoretical Determination of the Lift Coefficient 

for a Thin Elliptic Cylinder. Proc. Roy. Soc. (London), ser. 
A, voL 149, no. A168, Apr. 10, 1936, pp. 668-686. 

10. Goldstein, Sydney, ed.: Modern Developments in Fluid Dynamios. 

Clarendon Press (Oxford), 1938. 

11. Lamb, Horace: Hydrodynamics. Sixth ed., Dover Pub., 1932. 

12. Schlichting, Herman: Grenrsohioht-Theorie. Verlag und Druok 

G. Braun, Karlsruhe, 1951. (Available in English translation 
as NACA TM 1217.) 

TABLE I.— RESULTS AT SEPARATION POINTS 


Quantity 

Separation point (*=—0.1667) 

Top (7) 

Bottom (j) 

V 

80. 0° 

340. 83° 

R 

0. 985 

0. 364 

q 

1. 210 

-1. 204 

q. 

0. 139 

0. 137 


1. 192 

-7. 95 

<7» 

-0. 0062 

1. 49 

iT 

1. 11 

0. 41 

2 = 

13. 3 

-5. 6 

z. 

- 1. 405 

0. 52 

A' 

-55 6 

-5 4 

Y 

4 9 

-4 2 




